Interplay between Lattice Distortion and Spin- Orbit Coupling in Double Perovskites 
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We develop anisotropic pseudo-spin antiferromagnetic Heisenberg models for monoclinically distorted dou- 
ble perovskites. We focus on these A2BB'06 materials that have magnetic moments on the 4d or 5d transition 
metal B' ions, which form a face-centered cubic lattice. In these models, we consider local 2-axis distortion 
of B'-O octahedra, affecting relative occupancy of t2g orbitals, along with geometric effects of the monoclinic 
distortion and spin-orbit coupling. The resulting pseudo-spin- 1/2 models are solved in the saddle-point limit of 
the Sp(Ai') generalization of the Heisenberg model. The spin S in the SU(2) case generalizes as a parameter k 
controlling quantum fluctuation in the Sp(Af) case. We consider two different models that may be appropriate 
for these systems. In particular, using Heisenberg exchange parameters for La2LiMo06 from a spin-dimer cal- 
culation, we conclude that this pseudo-spin- 1/2 system may order, but will be very close to a disordered spin 
liquid state. 

PACS numbers: 75.10.Jm, 75.1().Kt 
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I. INTRODUCTION 

Geometrically frustrated magnets have been of great recent 
interest, and are a common starting point in search of exotic 
ground states.'-^ One class of such frustrated antiferromagnets 
is found in the double perovskite oxides, which host a wide 
range of interesting behavior ''"^ These compounds of chemi- 
cal formula A2BB'06 feature ordered, interpenetrating face- 
centered cubic (FCC) lattices of the B and B' ions when the 
charge difference between these ions is large. ^ Both B and B' 
transition metal ions are octahedrally coordinated by oxygen. 
A geometrically frustrated FCC lattice is obtained when only 
the B' ions are magnetic. 

A conventional picture of isotropic antiferromagnetic su- 
perexchange is insufficient for these materials. Altering this 
picture are two important effects considered in our work. The 
first effect is spin-orbit coupling, which is relevant for the 
Ad and bd transition metal ions that comprise the magnetic 
sites. Spin-orbit coupling has been seen to lead to increased 
correlation effects, particularly in materials containing hd Ir 
ions. This is responsible for topological insulating behavior,^" 
particularly in the pyrochlore iridates,""'* the Mott insula- 
tor ground state of Sr2lr04,''"^^ and the potential spin-Uquid 
ground state of Na4lr308^^"^^ and honeycomb compounds 
A2lr03.^'' Octahedral crystal fields favor the t2g c/-orbitals, 
which have an effective orbital angular momentum Lgff — 1, 
up to a sign difference. Combined with S = 1/2 spin angu- 
lar momentum, the pseudo-total angular momentum states of 
JeS = 1/2 and Jeff — 3/2 result. In this case, the quadru- 
plet of Jeff = 3/2 states form a lower energy manifold than 
the other two states of Jeff = 1/2. The second effect is ge- 
ometrical distortion from the cubic case; monoclinic distor- 
tion is commonly seen in double perovskites.^ Lowered sym- 
metry from the monoclinic distortion will spoil the exchange 
isotropy directly, and introduce new exchange pathways. One 
particularly important result is the local z-axis compression 
or expansion of the B'-O octahedra, which we refer to as a 
tetragonal distortion of these octahedra. While the octahedral 



crystal field favors the t2g orbitals over the eg ones, the tetrag- 
onal distortion will split the t2g levels. In the case of a local 
z-axis compression, the dxy orbital is favored to be occupied, 
while an expansion favors the d^z and dy^. All of these effects 
will generate the anisotropic interactions that form the focus 
of our models. 

The role of spin-orbit coupling in the undistorted cubic dou- 
ble perovskites has been carefully considered by Chen et al. 
for materials of d^ electronic configuration.^'* In this work, we 
focus on the Ad^ and 5c?^ monoclinically distorted double per- 
ovskites, and consider the quantum pseudo-spin- 1/2 models 
that result, as explained in the main body of the paper We are 
particularly interested in the case of a local B'-O z-axis com- 



FIG. 1: (color online). Magnetic ordering (Type I antiferromagnetic) 
of the spin-3/2 Ru in La2LiRu06 (blue, with arrows).'" Also shown 
are the non-magnetic Li (light green) atoms, and two of the Ru-O 
(purple) octahedra, showing the effects of monoclinic distortion. 
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pression, where orbital degeneracy is absent. La2LiMo06 is 
a candidate for such a material, while the otherwise isostruc- 
tural Sr2CaRe06 features instead a .z-axis expansion of the 
octahedra. La2LiMo06 shows no magnetic ordering down 
to 2 K from either heat capacity or neutron diffraction; how- 
ever, iiSR measurements show evidence of short-range corre- 
lations developing below 20 K.^^ The Curie- Weiss tempera- 
ture is negative, 9c = —45 K, indicating predominant anti- 
ferromagnetic superexchange. In contrast, Sr2CaRe06 shows 
spin-freezing behavior below 14 K.^' 

In the present work, we use the SpiN) generalization of 
Heisenberg models to describe these systems. This gen- 
eralization provides a unifying framework to study the effect 
of spin magnitude, from semiclassical ordering at "large spin" 
to possible spin liquid phases for "small spin". 

The ability to capture "large-spin" magnetic order may 
help to describe the higher-spin analogues of double 
perovskites. In particular, the "spin-3/2" analogue of 
La2LiMo06 is the isostructural La2LiRu06, whose con- 
figuration occupies all three t2g orbitals. Since the effective 
magnetic moment is close to the spin-3 / 2-only moment, there 
is only slight renormalization due to spin-orbit coupling, and 
intra-orbital Coulomb repulsion is the dominant effect in de- 
termining orbital occupancy. We model this material with a 
spin-3/2 Heisenberg model, given the lack of orbital degen- 
eracy, providing a test for Sp(A^)-predicted ordering at spin 
larger than 1/2. In fact, La2LiRu06 shows type I antiferro- 
magnetic ordering below 30 K,"^' where spins are aligned on 
each x-y plane but antiparallel on the x-z and y-z planes, as 
seen in Figure 1. This is consistent with the results in the 
semi-classical ("large spin") limit of our Sp(A^) model. In 
contrast, an appropriate pseudo-spin- 1/2 anisotropic Heisen- 
berg model for La2LiMo06 leads to the conclusion that this 
system must be very close to a spin liquid state. This may be 
consistent with the absence of magnetic order down to 2 K 
seen in experiment.'* 

The rest of the paper is organized as follows. In §11, we dis- 
cuss the effects of monoclinic distortion and spin-orbit cou- 
pling. This leads us to consider two different models, the pla- 
nar anisotropy and general anisotropy models, each taking the 
form of a pseudo-spin Heisenberg model. In §111, we solve for 
the classical spin ordering of both of these models. In §IV, we 
describe the Sp{N) generalization of the Heisenberg model 
and its mean-field treatment. Results of this mean-field treat- 
ment are shown in §V for the planar anisotropy model, and in 
§VI for the general anisotropy model. An extension to finite 
temperature is discussed in §VII. In §VIII, we summarize our 
results and discuss extensions of this work. 



II. MODEL 

In modelling monoclinically distorted double perovskites 
with Ad or 5d magnetic ions, there are two important effects 
of the monochnic distortion that should be considered in con- 
junction with spin-orbit coupling. The first effect of mono- 
clinic distortion is local z-axis compression or expansion of 
the B'-O octahedra, which affects orbital occupation. The sec- 



ond is the change of orbital orientation due to the geometric 
distortion, which affects overlap integrals and the resultant in- 
teractions. We will derive our models by considering the ef- 
fect of distortion and spin-orbit coupling on the interactions 
between t2g orbitals. 

One motivation for our models comes from a spin- 1/2 
Heisenberg model obtained via spin-dimer calculation for 
the isostructural monoclinically distorted double perovskites 
La2LiMo06 and Sr2CaRe06.'* In this method, the tetragonal 
compression (or expansion) of these materials was modelled 
by assuming occupation of only the d^y orbitals (or equal oc- 
cupation of only the d^z and dyz orbitals). This method is also 
sensitive to the effect of the geometric changes resulting from 
the distortion. However, spin-orbit coupling was not consid- 
ered, so that the assumed orbital occupation will be slightly 
incorrect. The result is an anisotropic 5=1/2 Heisenberg 
model, with estimates for the relative strengths of the cou- 
plings, seen in Table I. 



A. Interactions 

To understand the effects of the monoclinic distortion and 
spin-orbit coupling, we first look at the interactions between 
neighboring t2g orbitals in the case of cubic symmetry, as have 
been considered in detail by Chen et al?'^ To facilitate this, we 
show the six nearest-neighbor directions 6n for the FCC lat- 
tice in Figure 2. Without distortion, the a, b and c-axes are 
simply the Cartesian x, y and 2:-axes. The strongest interac- 
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FIG. 2: Nearest-neighbor lattice vectors dn and associated nearest- 
neighbor couplings J„ for the FCC lattice. 

tion is antiferromagnetic superexchange, involving sites and 
orbitals lying in the same plane. For instance, dxy orbitals on 
neighboring sites along the x-y plane will interact antiferro- 
magnetically. Ferromagnetic interactions between sites on a 
plane will couple orbitals lying on that plane to orbitals lying 
perpendicular to it."*"* Along the x-y plane, d^y orbitals inter- 
act ferromagnetically with neighboring dyz and dxz orbitals. 
Quadrupole-quadrupole interactions also exist between all t2g 
orbitals on neighboring sites, due to different orientations of 
the quadrupole moments of these orbitals. 



3 



Material 


Jl 


J2 


Ja 


J4 


Js 


Je 


La2LiMo06 


0.14 


1.0 


0.014 


0.014 


0.00043 


0.00043 


Sr2CaRe06 


0.87 


1.0 


0.16 


0.16 


0.25 


0.25 



TABLE I: Relative strengths of Heisenberg couplings, given in Fig- 
ure 2, from the spin-dimer calculation of Aharen et al?^ 



B. Monoclinic Distortion 

The first effect of monoclinic distortion is the local z- 
axis distortion of the B'-O octahedra, a compression for 
La2LiMo06 and an expansion for Sr2CaRe06. This splits the 
degeneracy of the three t2g orbitals. The d^z and dy^ orbitals 
will remain degenerate, but the d^y orbital will have a lower 
energy for a compression and a higher energy for an expan- 



sion. Consequently, the occupation of the d^y orbital will be 
favored or disfavored compared to occupation of the other two 
orbitals. This is taken as a very important effect in the spin- 
dimer calculation to explain the relative anisotropies of the 
two materials.^* 

The second important effect of the monocUnic distortion 
is a global c-axis elongation, and rotation of the B'-O octa- 
hedra, affecting the overlap of the occupied orbitals, which 
are now tilted out of plane. An example of this, in the case 
of La2LiRu06, is shown in Figure 1. The d^y orbitals, for 
instance, are tilted out of the a-b plane, and will have some 
interaction with d^y orbitals on neighboring planes. In this 
fashion, many new exchange pathways will contribute at the 
nearest-neighbor level. 

These effects generate a significant amount of exchange 
anisotropy in the spin-dimer calculation.^* The relative 
coupling strengths estimated by spin-dimer calculation for 
La2LiMo06 and Sr2CaRe06 can be seen in Table I. 
Interactions between x-y planes in La2LiMo06 are relatively 
weak, as expected from dominant in-plane dxy-d^y antifer- 
romagnetic interaction and c-axis elongation. We note that 
further in-plane anisotropy is significant, due to the strong ef- 
fect of Mo-0 octahedra rotation upon d^y orbital overlap. In 
Sr2CaRe06, intra-plane interactions are still larger than inter- 
plane interactions, even though the superexchange between 
dxy orbitals is not present. The only in-plane superexchange 
processes occur through tilted d^z or dyz orbitals. Neverthe- 
less, the inter-plane interactions are significantly stronger than 
in La2LiMo06. The length of the unit cell along the c axis is 
significantly larger than along the a or 6 axes, which could ex- 
plain the smaller inter-plane coupling compared to the intra- 
plane one. For both materials, however, the planar anisotropy 
of the couplings is clear, and effects of both geometrical dis- 
tortion and orbital occupation are important. 



C. Spin-Orbit Coupling 

Beyond monoclinic distortion, we now consider spin-orbit 
coupUng, which can be important in the 4rf and 5d magnetic 
ions commonly seen in the double perovskites. For instance, 
spin-orbit coupling in octahedrally coordinated Mo^+ is esti- 



mated to be on the order of 0. 1 eV.''^ The effect of spin-orbit 
coupling on the t2g orbitals of octahedrally coordinated ions 
is a well-studied problem. When the octahedral crystal field 
splitting is significantly large compared to the spin-orbit cou- 
pling, we may project out the Cg states. Upon projection, the 
L = 2 orbital angular momentum for the d orbitals looks like 
a L = 1 pseudo-angular momentum operator 1 up to a sign 
change, where L — >^ — 1. This Leff = 1 pseudo-orbital angular 
momentum combines with the S = 1/2 angular momentum 
of the single electron to create states of effective total angu- 
lar momentum Jeff = 3/2 and 1/2. The spin-orbit coupling 
AL • S breaks the degeneracy of these states, where the four 
Jeff — 3/2 states have an energy 3 A/ 2 lower than the two 
Jeff =1/2 ones. These J^s = 3/2 states are written in terms 
of the t2g ones as 
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With a d^ configuration, the occupancy of the dxy orbital upon 
projection to these states is given by^'' 



Tli, 



xy 



4 3^-^'^ 



(2) 



The occupation operators for the other t2g orbitals are given 
by cyclic permutation of the x,y,z indices, and the single- 
occupancy constraint ni,xy + rii^xz + ni^yz = 1 is satisfied. 

The effect of projection onto this Jeg = 3/2 subspace, due 
to large spin-orbit coupUng, has been considered by Chen et 
al. for the cubic materials.-'"* The Hamiltonian can be writ- 
ten in terms of the orbitally-resolved spin operators, such as 
^i,xy = SiUi^xy Upon projecting to the J^s = 3/2 states, 
these orbitally-resolved spin operators contain terms both lin- 
ear and cubic in j. The resulting Hamiltonian, containing 
terms of 4* and 6* order in j, leads to interesting multipo- 
lar behavior.-'"* 

When spin-orbit coupling is much larger than the local z- 
axis crystal field, the Jeff = 3/2 states provide the relevant 
starting point, rather than the t2g orbitals. However, one can 
consider the general splitting of t2g orbital degeneracy in the 
presence of both spin-orbit coupling and the local 0-axis dis- 
tortion. We can model each site with a local Hamiltonian 
-Hioc = A _ 2/3] - Al • S, where A > is the strength 
of the crystal field splitting due to local z-axis compression. 
The case for a local z-axis expansion has been considered 
by Jackeli and KhaliulUn.^^ We proceed in a similar manner, 
identifying the relevant low-energy eigenstates of Hioc- Di- 
agonalization of Hioc determines the lowest-energy Kramers 
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pair to be given by 



D. Planar-Anisotropy and General-Anisotropy Models 



|t)G = {i \yz, i) + \xz, i)) - I COS (Q) \xy, f) 

|4-)g = t) + Vz, X))-i cos {Q) \xy, i) 

tan(26») =2^A/(A + 2A). 



(3) 



The energy difference between the ground and first excited 
doublets is given by - A + ( A + 2 A) (1 + 1/ cos (26*) ) /4, which 
goes to zero as A 0, and approaches A— A/2 when A 3> A. 
We consider the case where this separation is large enough 
to focus on the lowest-energy doublet. This will require the 
tetragonal crystal field to be significantly larger than the ex- 
change coupling J, regardless of the relative strength of spin- 
orbit coupling. By projecting out the higher-energy states, we 
obtain a pseudo-spin- 1/2 model. 

Within this projection, we consider the form of the interac- 
tions in an otherwise cubic double perovskite, beginning with 
the quadrupole-quadrupole interaction. Due to the fixed or- 
bital occupation in (3), this interaction is constant and will not 
contribute to our models. The orbitally off-diagonal ferromag- 
netic interactions, of strength J', generate pseudo-spin inter- 
actions that are both spatially and spin-anisotropic. For our 
models, we will focus on the antiferromagnetic interactions. 
Nearest-neighbor interactions along the undistorted x-y, x-z 
and y-z planes are given by 

"Haf = J ^ f • Sj - 2 

<ij>inx-y 

+{xy yz) + {xy xz), (4) 

where rii^xy is the occupation operator of the dxy orbital at site 
i}'^ Upon projection to the lowest-energy doublet, we obtain 
a Heisenberg model in the pseudo-spin- 1/2 operators P,, 



<ij> inx-y 

^ sinW^^P,.P,+ ^ sinW^^P,.P,.. 



<ij> in x-z 



<ij> in y-z 



(5) 



For A ^ A, this result reduces to the one obtained by Chen 
et al. in the easy-plane limit of the cubic perovskite model 
with J' = 0.^"^ Without an accurate estimate for the strength 
of Hund's coupling to Coulomb repulsion, the ratio J'/ J is 
difficult to ascertain. However, we note that the easy-plane re- 
sult of Chen et al. is an antiferromagnetic state for J' < J?'^ 
Consequently, we consider the physical picture of antiferro- 
magnetic interactions, and as a first-order approximation we 
ignore the ferromagnetic contributions to the Hamiltonian. 

We note that the introduction of spin-orbit coupUng results 
in a reduction of the magnetic moment compared to the case 
of dxy occupation when A = 0. 



The first, and simpler, of the two models considered in this 
paper is concerned primarily with the effects of the tetrag- 
onal crystal field sphtting. Without spin-orbit couphng, we 
see easily from (4) that preferential dxy orbital occupation 
leads to anisotropic interactions that are stronger on the x-y 
planes. In this case, we have a true spin-1/2 antiferromagnetic 
Heisenberg model. However, considering spin-orbit coupling 
and tetragonal distortion leads to the pseudo-spin-1/2 antifer- 
romagnetic Heisenberg model in (5), with a similar form of 
anisotropy. From this, we are motivated to study the pseudo- 
spin- 1/2 antiferromagnetic Heisenberg model where coupling 
along the x-y plane differs from the couphng along the y-z 
and x-z planes. The planar anisotropy model is given in terms 
of pseudo-spin-1/2 operators (henceforth referred to as Si) by 



"HpA — Jin Si • Sj 

< ij> in x-y 

+ Jout Si • Sj + Jout Si • Sj. (6) 



<ij> in y-z 



<ij> in x-z 



Both Jin and Jout are antiferromagnetic, and one can consider 
this model as a generalization of the antiferromagnetic model 
in Eq. (5). The ratio Jout/ Jin depends on the strengths of the 
spin-orbit coupUng and tetragonal distortion of the octahedra, 
seen in A/A. In addition, it captures certain geometrical ef- 
fects of the monoclinic distortion, such as the global c-axis 
elongation, contributing to the particular planar anisotropy in 
(6). 

The other model considered in this paper will include in 
full the geometrical effects of the monoclinic distortion. This 
will generate many other anisotropic interactions, breaking 
the symmetry of the x-y plane. Effective pseudo-spin ex- 
change energies will become intrinsically anisotropic, in addi- 
tion to the effects of orbital occupation. We will model these 
hke the spin-dimer calculation does, with different strengths 
of the nearest-neighbor couplings shown in Fig. 2. Due to 
spin-orbit coupling, the particular parameters J„ in Table I 
will not be quantitatively correct. Nonetheless, we will con- 
sider them as a starting point to understand the effect of fur- 
ther anisotropy in the interactions. Estimates for corrections 
due to spin-orbit coupling are given in §VIB. The general 
anisotropy model is given by 



^GA = 5^ 5^ J„S(ri) • S(ri + 5n). 



(7) 



To analyze the model Hamiltonians (6) and (7), we will use 
the Sp(A'^) generalization of the Heisenberg model, which of- 
fers several advantages. The first is that the parameter TV al- 
lows for a controlled expansion, beginning from the saddle- 
point solution as A'^ — >^ 00. The second is that quantum fluctu- 
ations can be controlled by a parameter n (where n = 2S 
in the SU(2) case) allowing a transition from a classical- 
spin Umit (large k) to one dominated by quantum fluctuations 
(small k). This may capture a changing value of (pseudo)- 
spin. The gapped Z2 spin liquid, obtained as a disordered 



state in the Sp(A^) generalization, is often seen as a potential 
ground state in many Heisenberg models. ''■'•^'^ 

The Sp(A^) generalization may be capable of naturally 
capturing the changing behavior with S seen in the fam- 
ily of magnetic materials isostructural to La2LiMo06. The 
spin-3/2 La2LiRu06 is magnetically ordered, while spin-1/2 
La2LiMo06 shows short-range correlations and suppression 
of magnetic order. The isostructural spin-1 La2LiRe06 is 
more amenable to a multi-orbital model, and falls outside the 
scope of these calculations."*^ 



m. CLASSICAL ORDERING 



y 



^ — — — ^ 



In this section, we solve both planar anisotropy and general 
anisotropy models in the limit of classical spins. The mag- 
netic ordering patterns and wavevectors are determined by the 
OiN) model, where we generalize to A'' ^ oo components 
of the spin vector, as explained in Appendix A. We will see 
in §IV D that this corresponds also to the classical limit of the 
Sp(iV) model. 



A. Planar- Anisotropy Model 

In the planar anisotropy model (6), two phases are found 
with varying Jom/Jm, the ratio of inter-plane to intra-plane 
interactions. For Jout < Jin, the intra-plane interactions cre- 
ate antiferromagnetic Neel order within each x-y plane. For 
Jout > Jin, the inter-plane interactions create antiferromag- 
netic order between planes. 

For Jout > Jin, the ordering wavevector q is given by 



^(0,0,1) or ^(1,1,0). 



(8) 



Spins on each x-y plane are aligned, while spins on neigh- 
boring planes are antiparallel. Neel ordering is found along 
the x-z and y-z planes. The antiferromagnetic interactions 
between x-y layers are satisfied, as seen in Figure 3. 
For Jout < Jin, the ordering wavevector q is 

for arbitrary k^- Each x-y plane takes on the Neel order for 
a square lattice. The degeneracy in kz indicates that spins on 
neighboring planes may take any relative overall orientation. 
An example of this ordering, with k^ = 0, is given in Figure 
4. We will see in §IV D that this degeneracy is broken by the 
introduction of quantum fluctuations, choosing kz = 0. 

Both of these states show Type I antiferromagnetic ordering 
on the FCC lattice, where ordering is antiferromagnetic on 
two of the x-y, x-z or y-z planes, and ferromagnetic on the 
other. 



B. General Anisotropy Model 

The two parameter sets in Table I also yield antiferromag- 
netic ordering in the x-y plane, similar to the Jout < Jin case. 



FIG. 3: View along the 2-axis of FCC lattice magnetic ordering of 
the planar anisotropy model for Jout > Jin- The sohd lines indicate 
an x-y plane of the FCC lattice, while the dotted lines indicate a 
neighboring plane. Spins are aligned on each of the x-y planes, but 
Neel ordered along x-z or y-z planes. 



^ 



^ — ^ 



^ — ^ 



FIG. 4: View along the z-axis of FCC lattice magnetic ordering of 
the planar anisotropy model for Jout < Jin with kz = 0. The solid 
Unes indicate an x-y plane of the FCC lattice, while the dotted line 
indicates a neighboring plane. There is Neel ordering along each 
of the x-y and y-z planes, but ferromagnetic ordering along the x-z 
plane. Also possible is a state where the ferromagnetic ordering is 
along the y-z plane instead. 



However, the degeneracy of kz is broken here at the classical 
level, where kz = for both parameter sets. Ordering as in 
Figure 4 results. 



IV. SP(iV) MEAN FIELD THEORY 

A. Sp(iV) Generalization of the Spin Models 

The SpiN) method is a large-AT generalization of the 
Sch winger boson spin representation.^^" In the physical 
case N = 1, Sp(l) is isomorphic to SU(2), and we have 
the standard Schwinger boson representation wherein Sja = 

\^\ai'^a)afibif} and the boson number per site h\j3ia = ni, = 
2S determines the spin quantum number. Here, a, /3 =t, i la- 
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bel the primitive spin-1/2 species that comprise the full spin 
angular momentum. We generalize to 2N flavors of bosons, 
where a = (m, a), labelled by m = 1 . . . iV, and a J,, 
transforming under the group Sp(A^).^^ k = rib/N acts in 
analogous fashion to 2S in the SU(2) case, controlling the 
strength of quantum fluctuations. 

When generaUzed to Sp(A^), the Heisenberg Hamiltonian 
(7), up to constants involving ub, is written as 



(10) 



Here, J^i^ 
given by 



is a 2N X 2 A'^ block-diagonal antisymmetric tensor. 



Jr 



ma.m'cr' 




(11) 



B. Mean-Field States 

The quartic terms in (10) can be quadratically decoupled by 
the mean field 



Qin — 



N 



E 



(12) 



When the boson dispersion becomes gapless, we allow for 
a condensate, hn^ = Xia € C, where a =t, 4-. so that 
(biia) has a finite expectation value. This will account for the 
appearance of long-range magnetic order. 

The projective symmetric group analysis may be used to 
characterize possible mean-field ground states; for Sp(A^) this 
has been applied to many other Heisenberg models. Quali- 
tatively different states are distinguished by the value of a flux 
quantity for plaquettes of the lattice. The flux on a plaquette 
of sites a ... z is defined by the phase $ in''^ 



E' 



lab[ 



Ql,)...Qy.{-Qla). (13) 



C. Mean-Field Hamiltonian 

After decoupling in the site-independent Qn fields, the 
Hamiltonian (10) becomes 



H = Y,Jn 



N 



\m=2 



N 



N 

+h.c. + — \Qn\ 



E/^' + E ^IrHa^r + NX*,X^ 



7n=2 



(14) 



Here, the boson number constraint is enforced on average by 
the inclusion of the Lagrange multiplier /Zj. We assume tranl- 
sational invariance, with — ji. We have allowed the m = 1 
component to condense, represented by x1 € C. 

The saddle-point Hamiltonian (for N ^ oo) is de- 
rived in full in Appendix B. The first step is a Fourier 
transform defined by bi = -^/^ X^fe ^feC"*'^ ''*. The 
second step is a Bogoliubov transformation diagonalizing 
the Hamiltonian, yielding a quasiparticle energy ujk = 
\J ~ (X^n '^nQn sin {k ■ (5„))^. The transformation is de- 
fined by b = T~^7, where the Hamiltonian is diagonal in 
the 7 basis. The condensate enters only via the total density 



E 



, and ±ki, the wavevectors of the boson dis- 



persion minimum where the condensate forms. 
We then write the diagonaUzed Hamiltonian as 

^ = Eyl^?^l' + ^(-i-- + -) 

-Vn ^ Ji,Qi sin {k\ ■ 6) 

5 

k 

D. Semiclassical Lar^e-/t Limit 



(15) 



A nearest-neighbor Heisenberg model will favor the zero- 
flux states at small k, particularly for plaquettes of smaller 
length.'*^ On the bipartite cubic lattice, for instance, a transla- 
tionally invariant choice of Qij yields zero flux on any plaque- 
tte. Since the FCC lattice is frustrated, a translationally invari- 
ant Qij = —Qji, while giving zero flux on most plaquettes, 
leaves tt flux on a small number of plaquettes. In particular, 
assuming all Q to be translationally invariant and positive, the 
four-site plaquettes with tt flux have sites on both the x-y and 
y-z planes, such as i, i + ^i, i + i + ^3, where i and i ± y 
are joined by the plaquette. There are eight such plaquettes 
with TT flux, of a total of thirty-six four-site plaquettes involv- 
ing site i. This provides motivation to consider translationally 
invariant mean-field solutions, which we restrict ourselves to 
in this work. 



We take advantage of the Sp(A^) fluctuation parameter k to 
look at the semiclassical magnetic order from the k ^ 00 
limit. This provides a link from the classical order of §111 to 
the magnetic order seen at finite k. 

We begin by approximating the Hamiltonian for k » 1. 
Here, leading-order behavior in the Hamiltonian is of O(k^). 
Corrections, of 0(k), act to split degeneracy of the classical 
ordering. We have that Q, fi and n are all 0(k) as k ^ 1. 
Ec, the largest contribution to the energy is of 0(k;^): 



NsN ^ 2 



-l-n^ J5(55sin(A;i • (5), (16) 
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while the first-order quantum correction Ei, of 0(k;), is given 
by 
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El 



(17) 



where Q, fi and n are given by solutions minimizing the clas- 
sical energy (16).^^ The mean-field equations for Ec are eas- 
ily solved, yielding n = k, ^ = - J2n JnQn sin (ki ■ Sn), 
and Qm = — k sin (fci • Sm)- We can then write Ec as a func- 
tion of the minimum wavevector ki : 



Ec 



= _«2^:|isin2(ki.5„). (18) 



With the boson dispersion minimum at ±ki, spin ordering 
occurs at the wavevectors q = ±2ki. The minimum of Ec 
corresponds to an ordering pattern equivalent to that of the 
classical OiN) model (see Appendix A for details).''^ The cor- 
rection (17) can then easily be computed for all ki (with cor- 
responding Q, ji, n) in the degenerate set of minima of (18). 



V. PLANAR ANISOTROPY MODEL RESULTS 

In this section we study the planar anisotropy model with 
in-plane coupling Jin (Ji = J2) and out-of-plane coupling 
Jout ( Jz = J4 = J5 — Jq). We study the effect of quantum 
fluctuations, controlled by k, and coupUng anisotropy, con- 
trolled by Jout/Jin- In §in, we saw classical Neel ordering 
on each x-y plane. The first-order quantum correction Ei in 
(17) breaks the degeneracy. After this "order by disorder", the 
ordering wavevectors are 

Spins are aligned along either the x-z or y-z planes. Ordering 
along one such direction was seen in Figure 4. 

As K is reduced from this limit, we wish to see the evolu- 
tion of the ordering wavevector and mean-field parameters. 
For small k, we investigate the destruction of the ordered 
state by quantum fluctuations. We note that the semiclassi- 
cal solutions, for all values of Jont/Jm, all feature |Qi| = 
I (32 1, 1 I = IQ4I, and IQ5I = IQel- Motivated addition- 
ally by the equality of in-plane couplings, Ji = J2, and of 
between-plane couplings, J4 = J4 = J5 = Jg, we take an 
ansatz with Qi = Q2, Q3 = Qa, and Q5 = Qq. The relative 
signs, such as between Qi and Q2, correspond to making a 
particular gauge choice. With such an ansatz, the semiclassi- 
cal solutions remain unchanged, with wavevectors (8) or (19) 
as appropriate. Furthermore, relaxing the ansatz suggests that 
the equivalence \Qi\ = IO2I, IQ3I = IQ4I, and \Qr^\ = [Qgl 
is retained down to low k. With this ansatz, we numerically 
solve the mean-field equations, given explicitly in Appendix 
B. The resulting phase diagram is given in Fig. 5, in which 
there are five phases to consider. 
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FIG. 5: Heuristic phase diagram for the Qi,Q3,Q5 ansatz of the 
planar anisotropy model. Note that the label Qm = indicates that 
Qm is negligibly small (compared to k and the finite Q) in the con- 
densed phase; Qm is identically zero in the corresponding spin liquid 
phases. Solid lines indicate second-order transitions, while dashed 
lines indicate first-order transitions. 



A. Inter-plane Antiferromagnetic Order 

This state is an extension of the classically ordered state for 
Jout > Jin, with antiparallel magnetization on neighboring x- 
y planes. Ferromagnetic ordering is seen along the x-y plane, 
with Neel ordering along the x-z and y-z planes. In this state, 
the intra-plane Q 1 = <52 is significantly smaller than the intra- 
plane Q3 through Qe • The ordering wavevector has only smaU 
corrections to the classical result (8). 



B. x-y Plane Neel Order 

This state is an extension of the classically ordered state for 
Jout < Jin, with Neel order on the x-y planes. It is character- 
ized by large |Qi | = \Q2 \ within the x-y plane. Of the two in- 
dependent inter-plane Q, one is significantly smaller than the 
other, depending on the gauge choice of ferromagnetic order 
direction (along the x-z or y-z plane). The ordering wavevec- 
tor has only smaU corrections to the semiclassical result (19). 



C. Inter-Plane Spin Liquid 

This state is a disordered analogue of the inter-plane or- 
dered state (§VA) for Jout > Jin- However, the intra-plane 
Qi = Q2 are identically zero in this state. While the di- 
rect intra-plane correlations are consequently zero, the finite 
inter-plane Q prevent the lattice from decoupling. The mini- 
mum wavevector, determining short-range order, still has only 
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small corrections compared to the ordered minimum (8). The VI. GENERAL ANISOTROPY MODEL RESULTS 

transition into this state from the intra-plane ordered state, as 
K is lowered, is second-order. 

A. Spin Dimer Parameters 



D. Three-Dimensional Intra-Plane Spin Liquid 

This state is a disordered analogue of the inter-plane or- 
dered state (§VB) for Jout < Jin- However, one of the intra- 
plane Q is now identically zero, such as Q3 = Q4. The other 
intra-plane Q is nonzero, but still smaller than the in-plane 
Qi = Q2, preventing the lattice from decoupling. As be- 
fore, the minimum wavevector, determining short-range order, 
has only small corrections compared to the ordered minimum 
(19). The transition into this state from the intra-plane ordered 
state, as k is lowered, is second-order. 



E. Quasi-Two-Dimensional Spin Liquid 

In this state, all inter-plane Q vanish: = Q4 = 

Qb = Qe — 0. The system then consists of decoupled two- 
dimensional x-y planes in this mean-field theory. The tran- 
sitions into this state, from either the ordered or disordered 
intra-plane states for Jout < Jin, are weakly first-order. The 
minimum (short-range order) wavevector no longer takes the 
semiclassical value, instead taking a different value among the 
classical solutions (9), with kz ^ 0.15. 



F. IHcritical Point and Destruction of Order 

We find a tricritical point at Jout = 0.58 Jn separating the 
intra-plane spin-liquid phases from the x-y plane Neel ordered 
phase. For Jn > Jout > Jout, the ordered state first enters the 
three-dimensional spin-liquid state as k is decreased. A first- 
order transition to the two-dimensional spin liquid follows as 
K decreases further. The k range of this three-dimensional 
spin liquid narrows as Jout reaches tricritical point, as seen in 
Figure 5. For Jout < Jem, in-plane coupling pushes the system 
to decouple. However, we expect that the Q = decoupling 
seen in all three mean-field spin liquid states is an artifact of 
the mean-field theory, and that 1 /N corrections will restore a 
small yet non-zero value to these Q. 

The critical k value of the destruction of magnetic ordering, 
Kc, is fairly small in this planar anisotropy model. Kc ranges 
from 0.1 for large Jout to 0.4 for small Jout- In the physical 
N = 1 case, k = 1 corresponds to the "most quantum" limit 
of 5 = 1/2. Our N ^ 00 solution indicates that ordering is 
Ukely to occur, even though mean-field theory overestimates 
ordering. While Kc will differ in the exact N = 1 theory, 
the values of Kc ^ 0.1 — 0.4 are too small to account for the 
behavior of La2LiMo06. 



We now turn to the particular parameter set in Table 1 mod- 
elling La2LiMo06. We saw that the semi-classical limit led 
to Type I antiferromagnetic order, with Neel order on the x-y 
planes. As for the planar-anisotropy model, we take advan- 
tage of coupling symmetry to simplify the mean-field calcu- 
lation. We make the ansatz Qs = Q4 and = Qq, since 
J3 = J4 and J5 = Jg. The semiclassical result satisfies this, 
while relaxing the ansatz again suggests this structure carries 
to low K. Then we numerically solve the resulting mean- 
field equations. The mean-field solution finds that ordering 
persists down to Kc = 0.986. As in the planar anisotropy 
case, the ordering wavevector changes Uttle with k, and Q5 re- 
mains significantly smaller than the other Q. At Kc, there is a 
weakly first-order phase transition into a disordered state with 
Qi = Q3 = Q5 = 0. This highly anisotropic mean-field solu- 
tion consists of decoupled quasi-one-dimensional chains, with 
Q2 contributing the only non-zero correlation. The phase dia- 
gram for the general anisotropy model with parameters mod- 
elling La2LiMo06 is given in Figure 6. As before, we expect 
1/N corrections to remove this decoupUng. 



Anisotropic Spin Liquid AF Magnetic Order 
0.986 

FIG. 6: Phase diagram as a function of k for the general anisotropy 
model with parameters for La2LiMo06, from Table I. For n larger 
than Kc — 0.986, the system is in a three-dimensional magnetically 
ordered state, as in the semiclassical limit. For k smaller than Kc, the 
system is in an anisotropic and highly decoupled spin liquid state. 



The parameter set for Sr2CaRe06 in Table I behaves sim- 
ilarly, although the transition occurs at a smaller Kc — 0.41, 
similar to the values from the planar anisotropy model. 

Two comparisons to the planar anisotropy model are rele- 
vant. The first is that at large exchange anisotropy, the mean- 
field theory continues to predict immediate transitions from 
magnetic order into maximally decoupled spin liquid states. 
Additionally, this anisotropy stabilizes these decoupled states. 
For the La2LiMo06 parameters, we see a marked increase in 
Kc, which falls quite close to 1. This saddle-point solution 
suggests that the S = 1/2 system must be very close to the 
transition to a spin-liquid state, even if magnetic order even- 
tually appears at very low temperature. The effect of further 
quantum or thermal fluctuations may be sufficient to destroy 
the order. This could explain why no long-range order is ob- 
served in La2LiMo06 down to 2 K, while /xSR shows at most 
short-ranged order. The distortion of La2LiMo06 from the 
cubic perovskite structure is key in moving beyond the mag- 
netic order predicted by the planar anisotropy model. 
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B. Corrections to In-Plane and Out-of-Plane Anisotropy 

While the Table I parameters give a good picture of the 
anisotropy of La2LiMo06, they will not be quantitatively cor- 
rect. We wish to look at deviations due to the inclusion of 
spin-orbit coupling, from the viewpoint of in-plane and out- 
of-plane anistropy. The change in orbital occupation will re- 
sult in a reduction of (ij.j,-mediated coupling as spin-orbit cou- 
pling increases, along with new contributions, primarily out- 
of-plane, from dxz and dyz occupation. From these consid- 
erations, we estimate changes to J„ so as to minimize the re- 
sulting anisotropy, thus estimating a lower bound for Kc upon 
inclusion of spin-orbit coupling. We determine the effective 
couplings J„ in a manner similar to model (5), but with in- 
trinsically anisotropic exchange modified by orbital occupa- 
tion. In general, we have 

Jn ^ COS (efj^y + i sin {efj^'^y^, (20) 




Rl = 0.14 
Rl = 0.16 
Rl = 0.19 
Rj = 0.24 
Rj = 0.29 
R, = 0.34 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 

Ro 



FIG. 7: (color online). Critical value Kc of destruction of magnetic 
order. Two types of anisotropy are considered. Ri is the ratio of 
the anisotropy within x-y planes, while Ro is the ratio of anisotropy 
between these planes. 



with as defined in (3). While the 6 — spin-dimer parame- 
ters give J^y, the J^^'^^ are unknown. Since they arise from 
octahedral tilting, the in-plane J^^'^^ will be quite small, sim- 
ilar to how the out-of-plane J^^ are small. Since 0.25 sin [0)'^ 
is also small, we ignore that term by estimating Jf 2 = 0. 
For the out-of-plane interactions, we will make a large esti- 
mate for J^^^^y^ to minimize the out-of-plane anisotropy, by 
taking J3 4 l^g — J2 ^, the largest exchange scale in the prob- 
lem. In terms of the spin-dimer parameters J^^, we estimate 
the change in magnitude of Jn due to the change in orbital 
occupation from spin-orbit coupling by taking 



VII. FINITE TEMPERATURE 

Thermal fluctuations of the quasiparticles in (15) introduce, 
beyond quantum fluctuations, another mechanism inducing 
disorder At nonzero temperatures, these excitations have a 
thermal Bose distribution. The energy (H) and the mean-field 
equations, (15) and (B4), are modified accordingly. Thermal 
fluctuations will reduce magnetic ordering and correlations. 
We see different finite temperature behavior depending on the 
state (ordered or spin liquid) seen at T = for a given set of 
Jn and K. 



Ju2 = cos (9 fjl^, 

J3,4,5,6 = COS (0)^3^1^5,6 + ^ sin (0) ^3^^ . (21) 

For the case of A ^ A, we find that Kc reduces to 0.86. How- 
ever, for a moderate case of A = A, we find that there is only 
a slight reduction in Kc to 0.98. For moderate values of A/A, 
these mean-field results indicate that the system is still close to 
a disordered state; however, this will be sensitive to the value 
of A/A. 

Exchange anisotropy has shown to be very important, from 
the results for the spin-dimer parameters and the spin or- 
bit coupling rescaled values (21). To better understand the 
combined effect of in-plane and out-of-plane anisotropy, we 
consider a model with slightly less than the full anisotropy, 
where Ji = R1J2, J3 = J4 ~ R0J2, and J^ = Jq = 
RoRiJ2- This captures the in-plane (Ri) and out-of-plane 
(Ro) anisotropy, differing from the full anisotropy only in the 
very small exchange parameters J5 and Jg. In Figure 7 we 
show function of Rp, for several values of Rj. We see 

that Kc decreases fairly evenly as either Rq or Rj increases. 
This confirms that both in-plane and out-of-plane anisotropy 
are important in securing a large Kc- 



A. Zero-Temperature Disordered Phases 

From disordered phases, as T increases, the magnitudes of 
all Q decrease. The smaller the value of Q at T = 0, the lower 
the temperature at which Q reaches zero. At a large enough 
temperature, all Q are zero, describing a perfectly paramag- 
netic state, where spins are independent and completely un- 
corrected. This unphysical behavior at high temperature is 
typical of iV — > 00 solutions of Schwinger boson mean-field 
theories, and disappears for smaller values of N.*^ 



B. Zero-Temperature Magnetic Phases 

From ordered phases, as T increases, the condensate den- 
sity n decreases along with the mean-field parameters 
It similarly reaches zero at a large enough T. At large k, the 
transition to the perfect paramagnet state is first-order, with 
the system remaining in the ordered state until all Q and n 
discontinuously jump to zero. This occurs even for moderate 
values of k, such as k ^ 0.5 in the planar anisotropy model. 
For instance, with Jout = 0.54,/in and k = 0.5, this transition 
occurs at T 0.44 Jin. With Oc = -45 K and 5 = 1/2, 
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the transition temperature T = 53 K, an overestimate to be 
expected of mean-field theory. 

For smaller k, close to the disordered state boundary, the 
transition is second order. Furthermore, the order can be de- 
stroyed before the Q become zero; the system has a second- 
order transition to a thermally disordered state before entering 
the perfect paramagnet state. We show such an example in 
Figure 8. Here, k — 0.2, just above the zero-temperature crit- 
ical Kc for Jout = 0.54 Jin. At T = 0, the transition with vary- 
ing K went from ordered state directly into a quasi-two dimen- 
sional spin liquid. At finite temperature, we see that there is 



0.05 




-0.35 



0.05 0.1 0.15 0.2 0.25 0.3 

T(Ji„) 

FIG. 8: (color online). Mean field and condensate density (magenta) 
destruction with increasing temperature, shown for Jout = 0.54Jin 
and K = 0.2 in the planar anisotropy model. Above T ~ O.lJin, 
the magnetic ordering is destroyed, leaving a thermally disordered 
state. As Qs (blue) and Qi (red) become zero, the system enters a 
two-dimensionally or completely decoupled state, respectively. 



a window, O.lJin ^ T < 0.15Jin, where a three-dimensional 
disordered state exists, in contrast with the decoupling behav- 
ior of the r = mean-field theory. 

The general anisotropy model with La2LiMo06 parame- 
ters shows similar behavior. However, at k = 1, the transition 
from the ordered state looks weakly first-order, with the sys- 
tem directly entering a quasi-two-dimensional decoupled state 
where only Qi and Q2, both in the x-y plane, are nonzero. A 
fully three-dimensional disordered state is not predicted here 
by the finite-temperature mean-field theory. Nonetheless, this 
case illustrates how fluctuations destroy magnetic order and 
inhibit coupling in the spin-liquid states. As before, we ex- 
pect 1/iV corrections to further restore correlations. 



C. Heat Capacity 

The presence of the perfect paramagnet state is an artifact 
of the mean-field theory. Regardless, the magnetic contribu- 
tion to the heat capacity is an important physical quantity, and 
can be reliably calculated in this approach at low tempera- 
tures. Cv is found straightforwardly from d{H)/dT. In the 
magnetically ordered states, we find that Cy cx at low tem- 
peratures. This is expected from three-dimensional antiferro- 
magnetic spin wave contributions. In the disordered states, 
Cv oc exp (— Ac/fesr). Ag scales roughly with the spin 
gap, as expected for gapped states. Unfortunately, the lattice 
match material for La2LiMo6 was not useful in subtracting the 
lattice contribution to the heat capacity.-''' Without clear data 
for the magnetic contribution to the specific heat, direct com- 
parison is not feasible. For a system close to the ordering tran- 
sition, such as the general anisotropy model for La2LiMo06, 
the behavior persists only at extremely low temperatures, 
further complicating potential comparison. 



Vm. CONCLUSION 

We have modelled the effects of monoclinic distortion and 
spin-orbit coupling in 4d^ or 5d^ double perovskites. Local z- 
axis distortion of the magnetic ion-oxygen octahedra changed 
dxy orbital occupation compared to the other t2g orbitals. 
Geometrical effects of monoclinic distortion changed orbital 
overlaps, introduced multiple exchange pathways and gener- 
ated significant anisotropy. Considering spin-orbit coupling in 
conjunction with the local z-axis crystal field yielded a lowest- 
energy doublet of states and a pseudo-spin-1/2 Heisenberg 
model from antiferromagnetic interactions. We considered 
first the general case where interactions between sites on x- 
y planes differ in strength from interactions between these 
planes. This planar anisotropy model was studied for a gen- 
eral ratio of these two couplings. Geometrical changes of the 
monoclinic distortion induce further anisotropy among the in- 
teractions, especially within the x-y plane, leading to the gen- 
eral anisotropy model, studied for particular parameters mod- 
elling La2LiMo06, estimated from spin-dimer calculation.^^ 
We solved both these models in the saddle-point limit of 
the Sp(Af) generalization of the Heisenberg model. Semi- 
classical ordering was determined to be Type 1 antiferromag- 
netic, with antiferromagnetic order on two of the x-y, x-z, 
y-z planes, and ferromagnetic order on the other The Sp(7V) 
method connected the semiclassical results to the limit of large 
quantum fluctuations. The large interaction anisotropy of the 
general anisotropy model predicted disordering at a relatively 
large Kc — 0.986. The iV = 1 pseudo-spin 1/2 system was 
determined to be very close to a disordered state, even if or- 
der sets in at a low temperature. This could explain the lack 
of long-range order seen down to 2 K in La2LiMo06. Fur- 
thermore, estimates of the effect of spin-orbit coupling on the 
spin-dimer calculation parameters of Table 1 reduced Kc only 
to 0.98 for moderate strength of spin-orbit coupling. The sys- 
tem is still close to a disordered state in this case. 

Further experimental and theoretical inquiries follow as 
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natural extensions of our investigation. Single-crystal exper- 
imental results would be useful, primarily in determining the 
short-range ordering wavevector of La2LiMo06. Results at 
temperatures lower than 2 K could determine specifically how 
antiferromagnetic order is being suppressed. Finally, esti- 
mates of the strength of the spin-orbit coupling and crystal 
field splitting would guide a more precise model of the mon- 
oclinic distortion. 
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Appendix A: Classical 0(N) Model 

We begin by writing the real- space partition function for the 
Heisenberg Hamiltonian on the FCC lattice, 

Z = y £)0£)/iexp(— 5(0, /z)), where S'(0,/i) = 



The corresponding saddle-point solution gives from 



(Al) 



Here, the 0{N) model generalizes the spin 4> from a three- 
component vector to an N-component vector. The first 
step is to take the Fourier transform defined by = 
J2k "^k 6xp (— ik • Ti), where Ng is the number of sites 
of the lattice. After the Fourier transform, we have 




2 



N, 



o ■ E"^*"^°^('*'^") 

. n=l ) 



(A2) 
(A3) 



We perform the Gaussian integral over (/)k and i/)^, giving 



Z (X 



D 



D{k, m) = + E ' 

(A4) 



i = J-y I 



The spin-spin correlation function scales as 



(A5) 



(A6) 



As ,8 oc, the minimum of £)(k, /z) will become the 
dominant contribution; magnetic ordering will occur with the 
wavevector q that minimizes Xln=i "^n cos (q • (5„). 



Appendix B: Saddle-Point Solution 

To find the saddle-point solution, we first look at the Fourier 
transform, defined as bi = J^k bke~'^'''^\ After taking 
this transform, the Hamiltonian (14) becomes 



n 



E r + M f-i - « + ^ E ^^^^^1 

n \ ^ k J 

"aT E ( ^ —^aa'X%x'L,,e'^'' ^ + h.C. 



kn 



IN E (^l 



-Bk mA^-W' 
Bk=i^ JsnQsn sin {k ■ 5„) 

n 

(Bl) 

where Ng is the number of sites in the system. 

The quadratic part of the mean-field Hamilto- 
nian in (Bl) is diagonalized by a standard Bogoli- 
ubov tran sformation.^^ With the qu asiparticle energy 
ujk = [i^ - JnQn sin (k ■ (5„))2, the diagonaUzed 
quadratic terms are 

^ (l + 7i^7fcT + llpk^ ■ (B2) 
* k 

Here, the transformation is defined by b = T~^7, where 
the columns of are the eigenvectors of r/M, M is the 
quadratic Hamiltonian matrix in (Bl), and the 2N x 2N r/ is 
given by 



^a0 a < N 
— Sal3 a> N 



The structure of the condensate can be determined from the 
associated mean-field equation: d{H)/dx1 = 0. The solu- 
tion to the disordered case {x — 0) has a gapped dispersion. 
We can track when the gap vanishes and bosons begin to con- 
dense. We find that .t]^ is a linear combination of condensates 
at the minimum wavevectors ±ki: x\ = cidk-ki + C2Sk+ki- 
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We then rewrite the part of the mean-field energy depending In the condensed phase, to ensure a gapless dispersion, /i 



on and obtain the mean-field equation 



1 

n: 



JsQs sin {k ■ 6) 



— J2n JnQn sin (ki • dn) > 0. The form of x\. follows as 
x\. = -i(^5k-ki +ic*ih+ki- 



We arrive at the diagonalized Hamiltonian (15). From this 
follow the mean-field equations 



1 dE 

1 dE 

NsN dQm 
1 dE 



= -1- K+n+ —, 

N.^ujk 

= JmQm + TiJa sin (fci • 6m) 

= + ^ JnQn sin (fci • 6n) (if n > 0). 



1 ^ J2n JnQn sin {k- Sn) 

Ns ^ ujk 



{JmSin {k ■ Sm)) ■ 



(B4) 
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